/(/5 4 7002^ 


OCEANIC LITHOSPHERE AND AS THEN OS P HE RE : 
THERI-IAL AND MECHANICAL STRUCTURE 


f° 

'] 

I CD 

) 

j 03 

t 

CT> 

01 

«N 

«} cn, 

1 

H cnI 

M3 

V r-' 


C 3 cn. 

K! 

35 0? 


00 


■S3" 


Nf 


cnl 




f:3 t 

0 00 ( 

0 0 1 

» 0 


55 u m M 

H 

u . 

"a 



w a? tc 
w o 
W pq 

S 3 

m fT) 


0 

c> cn 

s 


EH 


H 


►-3 

Kl . 
«i > 

U 

S -H 

H 

03 c? 

ts; 

M ED 

rtt 

33 

W 

EH nJ 


•H 

0 

PC) P 5 
32 M 

Eh 0 



ir> 

'H 

cn 

•* H 

cn 

pq (tj 

CO 

eq U 


3; 
(^1 
J « CO 

o o 




CO 33 
! «; H 
CO 
’ «j 


P4 

« 

o 

U 

o 

« 

CO 


G. Schuibert 

Department of Geophysics and Space Physics 
University of California 
Los Angeles X California 90024 


C . Eroidevanx 

Labor atoire de Physique des Sol ides* 
Universite de Paris - Sud 
Orsay, France 91405 


D. A. Yuen 

Department of Geophysics an‘d Space Physics 
University of California 

V 

Los Angeles, California 90024 


September 24, 1975 
Revised 

February 4, 1976 








fee 












^Laboratoire associe an C.N.R.S 



2 


ABSTRACT 

We present a coupled thermal and mechanical solid 
state model of the oceanic lithosphere and asthenosphere 
which includes vertical conduction of heat with a temperature- 
dependent thermal conductiviry k(T), horizontal and verti- 
cal advection of heat, viscous dissipation or shear heating , 
and linear or nonlinear deformation mechanisms with temper- 
ature and pressure-dependent constitutive relations between 
shear stress and strain rate. A priori assumptions and 
inputs to the model are relatively few. We require a con- 
stant horizontal velocity u^ and temperature at the 
surface and zero horizontal velocity and' constant temper- 
ature at great depth. In addition to 'numerical values 
of the thermal and mechanical properties of the medium, 
v?e specify only the values of u^, and The 

model determines the depth and age-dependent temperature, 
hoxxzontal and vertical velocity, and viscosity struc- 
tures of the lithosphere and asthenosphere. In particu- 
lar, we deduce therefrom ocean floor topography, oceanic 
heat flow, and lithosphere thickness as functions of the 
age of the ocean floor. The model also yields the age- 
dependent shear stress in the lithosphere and asthenosphere. 
From the age-dependent geotherms and assumed values of the 
elastic parameters we construct seismic velocity profiles 
which exhibit a marked low velocity zone. To a first approx- 
imation, simple boundary layer cooling, i.e. the growth of 
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a cold thermal boundary layer by upward conduction of heat 
with k(T) and horizontal advection of heat with constant 
velocity u^, determines the thermal structure and quanti- 
ties derivable essentially solely therefrom^ namely oceanic 
heat flow, ocean floor topography, and seismic velocities. 
The simple boundary layer cooling solution gives heat flow 
«(age) ^ and depth of the ocean °=(age) for example. The 
second approximation, i.e. the self-consistent thermo- 
mechanical solution,, provides modifications to the thermal 
structure associated with shear heating, effects. This 
solution also gives the horizontal and vertical velocity, 
and viscosity structures not determined by simple boundary 
layer cooling and such quantities as the age— dependent 
lithosphere and asthenosphere thicknesses derivable from 
the mechanical field. Effects of viscous dissipation at 
old ages tend ro flatten heat flow, ocean depth and lith- 
osphere thicknesses, age curves, e.g. although ocean depth 
and lithosphere thickness tend to increase with age at 
all ages, shear heating tends to reduce the rate of increase 
at old ages below that given by (age) The importance 
of shear heatring increases with plate velocity, decreases 
with subasthenospheric temperature T^, increases with the 
activation volume describing the pressure-dependence of 
viscosity, and increases with age of the ocean floor. 
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INTRODUCTION 

The history of the ocean floor is now well documented 
from magnetic anomaly and deep sea drilling studies. It is 
generally acknowledged that geophysical characteristics like 
topography and heat flov/ are essentially a funcrion of the 
age of the ocean floor and/ therefore, are practically 
identical for a slow spreading ocean, the Atlantic, and a 
fast spreading ocean, the Pacific (Sclater and Francheteau, 

1970). This universal dependence on age is particularly 

striking for topography, whereas heat flow data are now 

known to be affected by hydrothermal circulation in the 

young oceanic crust so that the scatter in heat flow values 

is large and the average values of heat flow by conduction 

do not represent the entire heat flow to the surface (Lister , 1972 , 

1974) . For the Pacific ocean, the depths of the seismic 

low velocity channel for various ages are also available 

(Kausel et al., 1974). 

Physical models of the oceanic upper mantle providing 
a basic understrand-ing, of the geophysical data are derived 
from either the boundary layer cooling of a convective cell 
(Turcotxe and Oxburgh, 1967) or from the cooling of a rigid 
slab ox fixed thickness moving at constant velocity away 
from a hot vertical boundary (McKenzie, 1967) . The simple 
boundary layer cooling model leads to a /age proportionality 
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for quantities like depth of the ocean floor, inverse of 
surface heat flow, or depth to a given isotherm. For 
young ages, these proportionalities are also predicted by 
the rigid slab cooling model. Davis and Lister (1974) have 
in fact shown that ocean floor topography data for 0-80 My 
plot on a straight line vs. /age. For ages larger than 
its themai conduction time the constant thickness slab 
model of course shows departures from this simple behavior 
(Parsons and Sclater, 1976) . These authors have presented 
new topographic data corresponding to zero free air gravity 
anomaiies for the North Atlantic and for the North Pacific 
extending to 100 My and 145 My, respectively. The topographic 

data depart from a simple /age law for the oldest parts of' 

! 

the ocean floor. There is also a suggestion in the heat 

flow data of a flattening in the heat flux at great age 
(Sdsuar et al. , 1976). 

1 

These purely thermal models give, of course, no 
in mat ion about the mechanical structure of the upper 
manule, i.e. about the geometry of the lithosphere and 
asTiheno sphere and about the velocity field. A lithospheric 
thickness can however be defined if the assumption is made 
that a given isotherm, say 1200 °C, marks the transition 
between a solid lithosphere and an asthenosphere where 
partial melting generally occurs. This approach is based 
on a widespread belief that the seismic low velocity zone, 
which is well developed under oceanic regions, but not 
under old continental shields (Kngpoff , 1972) , corresponds 
to the occurrence of temperatures exceeding the solidus of 
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mantle material (Anderson and Sammis, 1970). 

In this spirit Parker and Oldenburg (1973) have cal- 
culated a lithospheric thickness proportional to /age. 

It should be stressed that this definition of the lith- 
osphere-asthenosphere boundary is different from that of 
plate tectonics which defines the asthenosphere as a zone 
of mechanical decoupling of the plates from the lower mantle. 
Furthermore / mechanisms other than partial melting may 
generally account for the zone of low velocity and high 
attenuation of seismic waves (Gueguen and Mercier, 1972) . 

It is however often assumed that this zone of low velocity 

corresponds roughly with the asthenosphere in the mechan- 

! 

ical sense. Physically, this coincidence could be correct 

I 

sines high temperatures can favor both dissipative phenomena 
(tfave attenuation, shear deformation) and decrease of elastic 
modn~* i. (wave velocities) . It is however not required that 

f 

parri al melting occur to have viscosities low enough to 
accoTint for lithospheric decoupling. Ample laboratory 
darca on the plastic deformation of olivine at high tem- 
perature and pressure have confirmed this (Post, 1973; 

Kohls tedt and Goetze, 1974) . 

In this paper we investigate both the thermal and 
mechanical structures of the oceanic upper mantle. Thermal 
quantities and associated topographic effects are computed 
simultaneously with mechanical quantities such as velocity 
and shear stress. Non-Nev/tonian and Newtonian rheological 
laws are used- We shall see that we account for available 
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geophysical data, in particular, for departures from a 
simple proportionality with /age. EKperimental values are 
used for the material properties so that in principle the 
model contains no us table parameters. Our results should 

help to formulate new geophysical programs aimed at a 
better understanding of the dynamic structure of tectonic 
plates . 
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QUALITATIVE DESCRIPTION OF THE MODEL 

The thermal and mechanical structures of the oceanic 
lithosphere and asthenosphere are coupled via the contri- 
bution of shear heating to the temperature distribution 
and the dependence of viscosity on temperature. Aside 
from the work of Schubert and Turcotte (1972) ^ we are un- 
aware of investigators who have considered this coupled 
thermo-mechanical problem; the point of view in most model- 
ling efforts has been to consider the material beneath the 
ocean floor to move horizontally with the constant velocity 
u^ and to solve the temperature equation by including hor- 
izontal advection of heat and upward conduction of heat 
(Langseth et al. , 1966; McKenzie, 1967; Turcotte and Oxburgh, 
1967) . In this paper we develop a theory for the coupled 
thermo-mechanical characteristics of the oceanic lithosphere 
and asthenosphere which includes important phenomena here- 
tofore largely" neglected in the modelling of these suboceanic 
structures ^ 

The model we consider is sketched in Figure 1. Ma- 
terial rising- from beneath the ridge moves essentially hor- 
izontaXIy away from the ridge, cooling to the surface as it 
ages. The cooling results in the formation of a boundary 
layer or lithosphere which is essentially rigid and moves 
with the horizontal velocity of -the surface u^. Beneath 
the lithosphere is the asthenosphere, a region wherein 
the mantle material is sheared and the horizontal velocity 
is reduced to zero. Within the asttencsphere there is also 
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a small upward velocity v to supply the mass required by the thicken- 
ing, horizontally moving lithosphere. The upward asthenospheric flow 
must be supplied by a deeper mantle circulation not considered here. 

On the left side of Figure 1 we list the quantities 
we input to our model. At the surface, we specify the 
temperature T = T^, the horizontal velocity u = u^ and the 
vertical velocity v == 0, At great depth we require the 
temperature to rise to 'a constant value and the hori- 
zontal velocity to be zero. The T^ boundary condition 
facilitates comparison with the simple boundary layer cool- 
ing solution. The right hand side of the figure shows 
the information we extract from the model calculations . 

This includes temperature, horizontal velocity, vertical 

velocity', shear stress, viscosity, and seismic velocities 

/ 

as functions of depth y and distance from the ridge x (or 

age — • From the thermal structure of the suboceanic 

manffce we also infer the topography , of the ocean floor and 

) 

surface heat flow as a function of age. 

Our model includes the following features ; 

1. Conduction of heat in the vertical, 

2. Thermal conductivity that depends on temperature. 

3. Advection of heat by a horizontal and vertical mass flow 

that depends on depth and age; the depth and age 
dependent vertical flow supplies the lithospheric 

mass. 

4. Contribution of viscous dissipation or ^ear heating 
to the thermal state, 

5. Linear or nonlinear deformation law relating shear 
stress and strain rate, 

6. Temperature and pressure dependent viscosity. 
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We derive the structure of the oceanic lithosphere 
and asthenosphere simply on the basis of subsolidus effects. 
This point of view contrasts with that of Parker and 
Oldenburg (1973) and Oldenburg (1975 ) , who define the base 
of the lithosphere as a solid-liquid phase boundary. 

As in our model of the continental asthenosphere 
(Froidevaux and Schubert, 1975) the shear stress will be 
found to take a nonzero value at the surface of the plate 
(y = 0) . This artificial stress at the surface replaces 
the driving forces acting, for example, at the plate boun- 
daries. Elsffiser (19 6^ has pointed out that the lithosphere 
may act as a stress guide. Our surface shear stress repre- 
sernrs this phenomenon. Its magnitude is not pre-determined. 

i 

Whar is imposed is the surface velocity. 
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MATHEMATICAL FORMULATION 

i 

As described in the preceding section, the suboceanic litho- 
sphere can be viewed as the relatively cold, near surface ther- 
mal boundary layer of a convecting mantle. As material spreads 
from an oceanic ridge, cooling to the surface occurs, resulting 
in the formation of the thermal boundary layer or lithosphere; 
the boundary layer thickens with increased distance from the 
spreading center, i.e.' with age of the lithosphere. The sit- 
uation has been sketched in Figure 1. 


Constant Thermal Conductivity and Horizontal Velocity 

The simplest mathematical description of the temperature 
profile in the oceanic upper mantle away from the ridge axis ex- 
presses a balance between conduction of heat vertically and advec- 
tion of heat horizontally, at a constant lithospheric velocity u^ 

whe‘re p is the lithospheric density, c is the specific heat at 

^ XT 

constant pressure, k is the thermal conductivity, x is the hori- 
zontal ' distance from the ridge, y is the depth and T is the tem- 
perature-. In writing the boundary layer equation (1) we have 

fi' ' ^ 

assumed! '<< ~ (i.e. vertical heat conduction is much more im- 
ox oy 

portant than horizontal conduction of heat) and p, c , u and k 

P o 

are constants. This boundary layer approach is not valid very close 
to "the ridge. It is an accurate description as long as the combined 
thickness of the lithosphere and asthenosphere is small compared with 
the distance from the ridge. Equation (1) is solved subject 


to the conditions that at the surface. 
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y = 0 , T is equal to and at gr eath depth , y / T 
tends to (a constant) . 

The solution of (1) , obtained by Turcotte and Oxburgh 
(1972) , is 

pu c , 

T = T^ + (T^-T^)erf{J(-~^) ^ , (2) 

kx 

where erf is the error function. The surface heat flux q 

^o 

(positive upwards) associated with this temperature dis- 
tribution is 

irkx 

As shown in (2) , the temperature depends on x and y in 

the combination y/vST. Surface heat flow decreases as 
-J' —3k 

X ^ or ^ge) and the "thickness" of the lithosphere grows 
as >6r or (age ) ^ , if lithospheric thickness can be defined 
as the depth to a particular isotherm. 

As the relatively cold lithosphere thickens with age, 
it contracts, and the weight of a vertical column of rock 
of constant height also increases with age. If we require 
isostasy, then the thermal contraction of the aging lith- 
osphere is reflected in the bathymetry, i.e. the depth of 
the ocean floor must increase with age. By accounting for 
the bathymetry of the ocean floor and the thermal contraction 
of the rocks with age, one finds equal masses in vertical 
col\3inns of rock and seawater if 
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D 


D 


ref 


a 



■/ dy(T-T^^£) 
o 


(4) 


where D is the depth of the ocean floor, a is the thermal 
expansivity, is rhe density of seawater and the sub- 

script ref indicates a reference condition, i.e. 

*^ref depth and temperature at a particular age of 

the ocean floor. In evaluating (4) , one may assume p is 
a constant. With the temperature profile (2 ) , we can inte- 
grate (4) analytically to obtain 



where the thermal diffusivity, is k/(pc^). Thus 
the depth of the ocean floor also increases as /x or v^age. 

This solution provides an adequate zeroth order des- 

r 

V 

cription of the temperature in the oceanic lithosphere. It 
does not, however, include the effects of a variable thermal 
conductivity , advection by a more realistic flow field 
and the viscous dissipation of the flow. Moreover, we have 
learned nothing about the structure of the velocity in the 
oceanic asthenosphere since we assumed, a priori , a con- 
stant horizontal velocity u^. In the following, we gen- 
eralize the solution (2)' for the temperature to account 
for the above-mentioned effects and simultaneously we 
determine the mechanical structure of the asthenosphere, 
wherein the horizontal velocity drops from its surface 


value to zero. 
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Variable Thermal Conductivxt:y , Constant Horizontal Velocity 
First we consider the effect of a thermal conductivity 
k which is an arbitrary function of temperature. In gen- 
eral, when k is not constant, equation (1) must be re- 
placed by 


3T 


d^T 3k 3T 
3y2 


( 6 ) 


In the case k = k{T) we can rewrxte.{6) as 


^ 3T _ , 3^T , ,3Tv 2 dk 
P’^o p 3x " ^ ^3y^ dT 


(7) 


This nonlinear equation is essentially no more difficult 

to rol-ve than (1) , the linear equation for constant k, since 

(7) also possesses a similarity solution. 

net k^ be the value of theimial conductivity at the 

su r-rac e and introduce the similarity variable 

\ 

5 = ‘ ( 8 ) 

'' k X 
o 


The dimensionless temperature 


0 


T-T 

o 



(9) 


then has to satisfy the ordinary differential equation 



(d?^ <d0 k^ 


( 10 ) 


and the boundary conditions 
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0 = 0 at 5 - 0 9 1 as 5 ->- «> . (11) 


The upward heat flux g, according to (8) and (9) , is given 
by 


q 


k (T -T^) / pu c" , 

3 '!' _ o «> o / o p k d9 

3y 2 V X ' k * d? 

o o 


( 12 ) 


and the topography of the ocean floor is calculated from 
(4) , (8) and (9) . It becomes 


D 


D 


ref 




/ (1-6 (^))d? 
o 


(13) 


Thus heat flux^ bathymetry, and lithosphere thickness depend 
on age in theisame manner as they do for constant k; only 
proportionality factors are modified by variable k . 

We have nume-rically integrated (10) and (11) for the 
experimental Thermal conductivity function of Schatz and 
Simmons (1972) 

s 


^ ^ 88.33 

k 31+0 .-21T 
o 


for T < 500 K 


and 


k. 


88.33 
31+0. 2 IT 


+ 0 


485x10 ^(T-500) 


for T > 500 K 


(14) 
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where T is in degrees Kelvin and is 0.0113 cal/ (cm sec K) 
for an assumed surface temperature of 0 °C. 

The solutions show that the conductivity decreases 
with 5 ? small (conductivity decreases with depth 

near the surface) and then increases gradually toward the 
value k(T ). For T = 1200 the minimum value of k/k 
is about 0.6 at ?~0.4 and the asymptotic value is 0.73. 

For = 1600 °C the minimum is also about 0.6 at ^ ~ 0-5 and the 
asymptotic value is '•0.88. Within -the upper mantle there 
is about a factor of 2 variation in the thermal conductivity - 

The surface heat flux from (10) - (12) and (14) is 


pu c , 

= c^k^(T„-T^) (-,^^)^ 




(15) 


where c. is 0.4656 for T = 1200 °C and 0.4680 for T = 

JL CO CO 

1600 °C. Since (15) is an exact solution- for variable k, 
it is not necessary to use (3) together with some estimate 
of an average k. to account for variable conductivity. 
Finally, we note that the bathymetry, from (10) , (11) , (13) 
and (14), yields 



where is 0.932 for T = 1200 and 0.948 for T = 

1600 °C. 

The temperature profiles derived from this similarity 
treatment will be shown as reference profiles later on 
(see e.g. Figure 4). 



17 


T^vto— D xiueiis ional Shear Flov?^ Viscous Dissipation./ Variable 
Thermal "Conduc tivi by ~~ 

Exact similarity solutions are no longer possible when 
one considers the combined thermo-mechanical structure of 
the oceanic lithosphere and asthenosphere. We formulate 
instead a quasi-one-dimensional or quasi-similarity approach 
to the problem which uses the simple boundary layer cooling* 
solutions previously discussed as a starting point. The 
model, as previously outlined in connection with Figure 1, 
includes the two-dimensional flow of a medium with a Newtonian 
or non-Newtonian temperature and pressure dependent rheology 
forced to deform with horizontal velocity u^ at the surface 
and xero horizontal velocity at great depth. The depth 
interval in which there is shear in the horizontal velocity 
u Is" the as theno sphere. The asthenosphere is simply a zone of 
decoupling for the plate; it is not 'required that a shallow return 
flow; occurs within its depth range. Conservation of mass requires 
a vertical transport with velocity v to balance the variation in 
horizontal mass flux due to the change in u(x) with age (x/u^) 
or distance from the ridge x. The temperature equation includes 
advection of heat by both the horizontal and vertical motions, the 
viscous dissipation of the flow, and variable thermal conductivity 

We are interested in describing the structure of the oceanic 
lithosphere and asthenosphere which extends over horizontal dis- 
tances of 1000 ~ 10,000 km but is only 100 to 300 km thick; thus 
we have again y/x << 1 and 3/3x << 3/3y. With the Boussinesg 
assumption, conservation of mass is 
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8u 3v 
3x 3y 


0 


(17) 


v/u is 0 (y/x) and u-^ is 0 (v|^) . These order of magni- 
tude relationships permit us to write the temperature 
equation in the approximate form 


, 3T 3Tv 

pc (utt- + 
p 3x 3y 



>3T. 2 3k 
^3y^ 3T 



(18) 


where t is the shear stress acting parallel to a horizon- 
tal surface element in the medium. The third term on 
the right hand side of (18) is the dominant contribution 
to shear stress heating or viscous dissipation. 

The rheological law connecting shear stress and strain 
rate takes the simplified form 


3~g 

5y 


T 


n 


T 


exp{- 


(E*+pV*) 

RT 


} 


(19) 


where. E* is the activation energy, V* is the activation 
volume, p is the pressure, R is the gas constant and is 
a -proportionality factor. For Newtonian flow n = 1 and 
xorr non-Newtonian flow n is larger than 1- We will assume 
in the following that n = 3 for non-Newtonian flow, a value appro 
priate for the deformation of olivine. The T~^ factor is not re- 
solved in experiments , but derives from theoretical models of dis 
location motions . 

If we assume infinite Prandtl ntimber so that the momentum 
equation expresses a balance between pressure, viscous and 
gravity forces and further assume zero horizontal pressure 
gradient, then we may approximate the pressure by its hydro- 
static value and the horizontal momentum equation becomes 
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9t 

3y 


0 


( 20 ) 


We take the horizontal pressure gradient to be zero since 
we are mainly interested in the relatively simple case of 
asthenospheric decoupling of a rigid moving lithosphere 
from a quasi-static deep mantle. A basically horizontal 
shear flow can accomplish this decoupling without the 
imposition of an external horizontal pressure gradient. 

We will consider the effect of a non-zero pressure gradi- 
ent in driving a return flow in the oceanic asthenosphere 
in a subsequent paper. The boundary conditions for (17) - 
(20) are 

Uq/ T = T^, V = 0 on y = 0, 

U 0',/T aS"Y ^ 00 (21) 


Equations (17) - (20) are very similar to the ones 
investigated by Froidevaux and Schubert (1975) in their 
study of the continental asthenosphere. The equations 
for the continental asthenosphere were strictly one- 
dimensional, i..e. V = 0 and = 0. However, the cooling 
of the oceanrc: lithosphere as it moves away from the ridge 
and ages' req uir es that we include the advective terms on 
the left hand side of the temperature equation in our des- 
cription of the oceanic asthenosphere. The equations are 
also similar to the ones used by Schubert and Turcotte 
(1972) in their one-dimensional analysis of a return flow 
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in the oceanic asthenosphere. Schubert and Turcotte (1972) 
did not include the advectxve terms 'in the temperature 
equation but they partly accounted for lithospheric cool- 
ing by integrating the one-dimensional equations with respect 
to depth, starting at an a priori specified depth where the 
oceanic temperature was assumed to already have risen to a 
sufficiently high value, over 800 *^C, that flow readily 
occurred, i.e. they began their integration essentially at 
the base of the lithosphere whose depth and temperature 
were assumed. In this paper, however,, we integrate the 
equations from the surface of the ocean floor downward 
and determine the thermal and mechanical structures of the 
oceanic lithosphere and asthenosphere self-consistently. 

We solve the system (17) - (21) by a method of successxve 
approximation which can, in principle, be carried through to 
an exact solution of the two-dimensional problem. In the 
first step, we assume that the form of the simple boundary 
layer cooling solution for the temperature structure of 
the lithosphere (henceforth when we refer to simple boundary 
layer cooling we m.ean the case u = u^, v = 0 and variable k) 
provides a wey ot relating local horizontal derivatives to 
local vertical derivatives. Since boundary layer cooling 
with u = u^ and v = 0 gives a thermal structure dependent 
on y and, x_, in the combination y//x", we initially take 


= _ i I ^ 

3x 2 X 3y 


( 22 ) 
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Then the system (17) - (21) reduces to a set of ordinary non- 
linear differential equations with y as the independent 
variable and x as a parameter 


y du _ ^ 
2 X dy dy 


(23) 


/ 1 y \dT 

pc_ (v - *^u) 

^ p 2 X dy 


= k- 


2 

d T 

^ 2 
dy 


+ 


, dk > , dT , 2 
^dT^ ^dy^ 


+ T§ , (24) 


du 

dy 


2B 


n n , 
exp(- 


(E*+pV*) . 
RT ’ 


(25) 


dX 

dv 


0 


(26)'- 


Because (23) - (26) are coupled equations with the 
parameter x/ their solution subject to conditions (21) will 

t ' 

depend on distance from the ridge, i.e. on age. The solution 
of (2J.) and (23) - (26) thus provides information on the age 
dependences of the thermal and mechanical properties of the 
lithosphere and asthenosphere not assumed a priori . A 
second approximation can further improve the solution. In- 
stead of replacing 3T/3 x and 3u/9 x in (17) - (21) by the 
approximation (22), we can evaluate these derivatives, as 
functions of x and y, from the solutions of (21) and (23) - 
(26) since we can obtain these solutions at arbitrarily 
close values of x. Equations (23) - (26) would then be 
replaced by a more accurate set of ordinary differential 
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equations in y. By repeating this process a number of 
times, a solution of the two-dimensional problem can be 
obtained to any desired accuracy. In this paper, we obtain 
the solution to the first approximation only, i.e. we solve 
(21) and (23) - (26). 

Equation (24) contains a second derivative of T; the 
two boundary conditions on T, namely T = on y = 0 and 
T -> T as Y ^ “ are thus sufficient to determine the tern- 
perature. Since (23) - (25) contain only first derivatives 
of velocity, the boundary conditions u = u^ and v = 0 on 
y - 0 suffice for the solution of the velocity field. From 
(26) we know that the shear stress is independent of depth. 
Howeveir, its value is not known a priori and further, it is 
a iuncuion of jt. The condition u 0 as y “ is an addi- 
tional. constraint on the solution of (21) and (23) - (26) 
that serves to determine t. Thus we not only find the 

X ♦ • 

t emper ature and velocity fields by splving this system, but 
in addition we solve for the a -priori unknown shear stress 
and its dependence on age. In connection with the age de- 
pendence of T , we note that if the pressure is to be approx- 
imated by the lithostatic pressure then 


dy 


0 (pg) 


>> 


dT 

dx 


where g is the acceleration of gravity. That this is an 
excellent approximation can be seen a posteriori from the 


solutions for x. 
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We will be interested in the profile of viscosity in 
the asthenosphere and the variation -of the viscosity struc- 
ture with age. Dynamic viscosity n is related to temper- 
ature and pressure by 


n = 


2e 


2B^t 

Di 


n-1 RT ' 


( 27 ) 


where e is the strain rate — ^ Note that n depends on 

T for n 1^ i.e. that any viscosity value refers’ to a given 
stress level or rate of deformation for non-Newtonian rheologies 
One speaks of an effective viscosity for such nonlinear rheolo- 
gies. 
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SOLUTIONS TO THE THERMO-MECHANICAL MODEL 

Temperature and velocity profiles together with the 
value of the shear stress have been calculated on the basis 
of -the set of equations and boundary conditions given above 
and for parameter values eniimerated below. This section 
describes the results both for non-Newtonian and Newtonian 
rheologies and emphasizes the main chacteristics of verti- 
cal structure of our model mantle at 2 ages, 10 My and 
150 My. The variation of the lithosphere-asthenosphere 
structure with age will be presented in the next section. 

Parameter Values 

Our numericaJL solutions are based on the following 
values of thermal, mechanical , and kinematic parameters. 

For the surface, temperature T^ we take O^C and for the 
temperature at great depth T^ we consider values between 
1200 *^C and IHOO °C. We have already noted that the value 
of thermal conductivity at the surface is 0.01132 cal/cm sec K, 
while (14) gives the value of k at temperature T. The 
pressure- dependence of k is not included. The specific heat 
at constant ■pr-essure c^ is assumed to be 0.27 cal/ (g K) and 

3 

the density p is 3.3 g/cm . Thus the value of the thermal 

diffusivrty at the surface temperature, dc~' 

—2 2 ^ 

1.27x10 cm /sec. The coefficient of theirmal expansivity 

a is somewhat uncertain and furtheirmore is temperature 

dependent. Its value enters pur analysis only in the 



25 


calculation of the topography of the ocean floor? we con- 

~5 “1 -5 ‘-r 

sider averaged values of a between 3x10 K and 4x10 K 
(Clark, 1966) . To encompass the range of observed plate 
velocities xn our calculations we take u^ equal to 2 cm/yr 
and 10 cm/yr, representative of the North Atlantic and 
Pacific oceans, respectively. 

The extent to which the simple boundary layer cooling 
solutions are modified by two-dimensional shear flow and 
viscous dissipation therein depends on the type of rheology 
and rheological parameter values assumed for the deforma- 
tion. We have investigated both Newtonian (n = 1) and 
non-Newtonian (n = 3) rheologies. In each case, two sets 
of xheological parameter values were considered, correspond- 

t 

xng ro materials which deform relatively easily and those 
whicni are relatively harder to shear. 

The rheological' parameter values are summarized in 
TabLe= 1. Since the activation volume V* is still an un- 
certain parameter we have considered it in the range 11- 

3 

30 cm /mole. It should be noted that the rheological 
parameter values for the Newtonian cases are based on 
quasi-theoretical estimates . Laboratory results for oli- 
vine lead to a nonlinear rheology with n = 3 . The olivine 
data of Table 1 were not all available when the continen- 
tal upper mantle model was investigated (Froidevaux and 
Schubert, 1975 (see, in particular, the note added in proof)). 
Thus the value of B used here for dry olivine differs from 
that used in our earlier paper. Another difference from 
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our continental model is the absence of both radioactive 
heat sources and heat flow from great depth. Such charac- 
teristics could be included in future extensions of the 
oceanic model. 

Temperature Profiles in a Dry Olivine Mantle 

We present the solutions of (23) - (26) with boundary- 

conditions (21) for a number of cases based on the dry 

olivine rheological parameter values given in coltamn 1 of 

Table 1. The solutions are for two values of u , T , and 

o oo 

V* and two different ages. 

Eirst we concentrate on the temperature vs. depth 
curves for the case T^ = 1600 '^C. The top of Figure 2 
shows two temperature profiles for a plate moving at 
2 ciri/yr- with respect to the lower mantle, one for an age 
of 10 3ly, the other for. an age of 150 My. These profiles 

I 

ane essentially identical to the simple boundary layer 
cooling solution (7 ) , so that viscous dissipation is negligi- 
ble^ This is still nearly the case when the plate velocity 
is 10 cm/yr, as illustrated at the bottom of Figure 2. 
Rowever, here the 150 My temperature curve, in fact, lies 
slightly above the corresponding boundary layer cooling 
solution (not shown on the figure ) , so that the heat flow 
at the surface is 1.13 HFU as compared to the boundary 
layer cooling result of 1.09 HFU (see Table 2). 
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Such dissipative effects are more pronounced xn Figure 3, 
which is simxlar to Figure 2, but has = 1200 °C, i.e. 
a cooler region at great depths For an age of 10 My the 
temperature profile is still identical to the boundary layer 
cooling solution for a slow plate (top of figure) but departs 
slightly from the bpundary layer cooling solution {dotted 
curve) for a fast plate (bottom of figure) . At 150 My shear 
heating represents an important contribution to the heat 
balance and the temperature profile departs significantly 
from that of the simple boundary layer cooling case (dotted 
curve) . . This introduces a dependence of the thermal struc- 
ture on the activation volume V* . The temperature profiles 
for V*’ between 11 and 3*0 lie in the shaded regions 

between the. curves appropriate to these limiting values of 
V*. Higher activation volumes imply higher temperatures for 
a given depth- Nocice 'finally that the bump on the T-curve 
is more pronounced for the fast plate (bottom of Figure 3) 
than for the slow plate (top of Figure 3) . At 150 My and 

T = 1200 ^Cf the calculated surface heat flow is thus 
00 

larger than the" value 0.82 HFU for boundary layer cooling. 
Depending on V* and plate velocity, can be as large as 
1.08 HHJ (see Table 2). 

We have thus seen that for certain plate ages and 
velocities, dt is physically possible to have a steady state 
temperature vs. depth curve which exhibits a change in the 
sign of its curvature. This, in a one-dimensional model. 
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would only be possible if heat sinks were present, if the 
thermal conductivity k had a peculiar T-dependence or, 
or course, if the situation were not steady state. Here, 
this phenomenon is related to the steady state advection 
of colder material from below in the two-dimensional flow 
field which will be described shortly. 

The importance of shear heating in the oceanic asthen- 
osphere is also illustrated by the fact that the functional 
dependence of on (T^ “ longer the linear one 

predicted by the simple boundary layer cooling solution 
(see (15) ) . This is shown in Figure 4 where the calculated 
value of at 1000 km from the ridge is shown vs. (T - T ) 
for plate velocities of 2, 5 and 10 cm/yr. The lower branches 
of the curves refer to the non-Newtonian cases discussed so 
far and the Newtonian I rheology (Table 1) ; the upper 
branches are far the Newtonian II rheology (Table 1) . Above 
about 1400 °C the linear relationship holds almost exactly, 
i.e. viscous dissipation is negligible. For larger dis- 
tances from the ridge we have seen that shear heating is 
noticeable even for T = 1600 °C (see q in Table 2) . Thus 
the nonlinearity would be much more noticeable for older 
ages than those of Figure 4, 

Depth of the Decoupling Shear Zone 

Horizontal velocity distributions u(y) corresponding 
to the above described temperature profiles for a dry 
olivine mantle are also shown in Figures 2 and 3 . The 
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associated viscosity profiles ri (y) are given in Figures 5 
and 6, respectively; calculated values of the shear stress 
T are listed in Table 3 . 

In Figures 2 and 3 we see that the horizontal velocity 
drops from its surface value of 2 or 10 cm/yr to zero in 
a narrow zone of shear, which we can also call the asthen- 
osphere since it corresponds to a region where the viscosity 
reaches a minimum as seen in Figures 5 and 6. The rather 
strong rise in viscosity below the asthenosphere is re- 
lated to the particular boundary condition we impose at 
great depth, i.e. T = T^. For a finite heat flow from 
great depth, 0 , the viscosity would rise less steeply 

belxaw the asthenosphere. 

The shear zone corresponding to a given geotherm is 

located in a depth range which depends upon the value of 

V", even when the geotherm is independent of V*. For 

1 

3 I 

V- — 30 cm /mole the reduction in horizontal velocxty takes 

3 

piece at shallower depth than for V* = 11 cm /mole. Since 
aTI solutions for intermediate values of V* fall in between 
these two extreme velocity profiles, the included area has 
been shaded to indicate the range of possible solutions due 
to the uncertainty in the value of V* . The' observed up- 
ward shift of the shear zone when V* increases, everything 
else being equal, simply indicates that when pressure effects 
in (25) or {27} become important, the system adjusts by 
moving up the shear zone. The depth of the shear zone is 
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a compromise between the opposing effects of temperature and 
pressure which tend to move the shear zone downward and upward, 
respectively. Thus in a given situation, the actual depth of 
the asthenosphere is determined by the particular values of 
the rheological parameters. 

The width of the shear zone can be defined as the depth dif- 
ference between u = 0.95 u and u = 0.05 u . From the figures we 
see that the widths range from 60 to 150 km for an age of 10 My. 

For a given value of the rheological parameter V*, the width 
of the asthenosphere does not depend very strongly on the plate 
velocity or age. This is apparent from the figures since the 
respective u profiles are almost parallel to each other. One should 
realize that the boundary condition g^ = 0 at great depth has an 

t 

influence on the location of the bottom of the shear zone. A nonzero 
value of is expected to extend the asthenosphere to greater depth. 

The depth of the top of the shear zone, i.e. the plate 

f 

thickness, does not vary significantly with plate velocity, 

\ 

as can be seen by comparing any pair of corresponding curves 
u(yi) for 2 and 10 cm/yr from the top and bottom parts of 
Figures 2 or 3. The plate thickness, on the other hand, 
exhibits a variation with age which will be analyzed in the 
next section. Here we notice, by inspection of Figures 5 
and 6 , that the value of the viscosity in the asthenosphere , 
i.e. in the shear zone, increases by up to one order of 
magnitude as the plate age increases from 10 to 150 My . For 
example, at a velocity of 2 cm/yr and V* = 30 cm /mole we 
find, from Figures 5 and 6, that the viscosity minima for 
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= 1600 °C are 5.5x10^^ and 2.6x10^^ poises, for 10 
and 150 My, respectively, and for = 1200 °C the 
viscosity mrnima are 1-2x10 and 1.2x10 poises for the 
two ages. This trend with age is less pronounced in cases 
where shear heating gives rise to strong departures of 
the old age geothe3rni from the simple boundary layer cool- 
. ing solution. This is seen for the fast plate with T^ = 

1200 °C at the bottom of Figure 6. 

The values of the shear stress t obtained for the 16 
cases discussed so far are listed in Table 3- The increase 
with age of the viscosity in the shear zone is seen to 
correspond' to an increase of the shear stress with age. 

For instance, in the two examples just discussed, one sees 
-that T varies from 0,90 to 17.0 bars in. the example of 
Figure 5 (T^ = 1600 ®C) and from 24 to 193 bars in the 
example of Figure 6 (1200 °C) , as age increases from 10 to 
150 My. The la?tter value of x is high and corresponds to 
a cold mantle (T = 1200 and an old (150 My)-, slow 

to 

3 ‘ ' 

(u^ = 2 cm/yr:) plate with a stiff rheology (V* = 30 cm /mole) 
Having 'this same plate move at 10 cm/yr generates more 
shear heating and x drops to 75.5 bars. One notices the 
relatively lower magnitudes of the shear stresses for 
T^ - 1600 °C as compared with T^ = 1200 °C. 

In situations where the entire heat production inside 
the model mantle is due to shear heating (Schubert and 
Turcotte, 1972; Froidevaux and Schubert, 1975) the shear 
stress X is inversely proportional to the plate velocity 
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u if- the heat flow is kept constant. Indeed, this heat 
o 

flow is equal to the energy dissipated tu^ in. a steady 
state regime. As heat is xntroduced by radioactive sources 
or as in this paper, by advection, this simple behavior is 
no longer true. Table 3 illustrates this point. An in- 
crease of velocity by a factor 5 , everything else being 
equal, produces an increase in t by a factor <2. When 
strong viscous heating effects are present (T^ = 1200 °C, 
age of 150 My) the shear stress decreases by a factor 1.4 
to 2.5. The assumption that t is proportional to u^, used 
by various authors (Solomon and Sleep, 1974; Forsyth and 
Uyeda, 1975; Solomon et al., 1975) in their discussions of 
plate motion and of the possible driving forces, is thus not gener- 
ally supported by thermo-mechanical models of plate de- 
coupling. 

Vertical Velocity Fields 

The vertical velocity v(y) , taken as positive upward, 

is also illustrated in Figures 5 and 6 for the 16 cases 

based on a dry olivine rheology. For an age of 10 My, v 

increases from z-ero just at the top of the asthenosphere 

and approaches an asymptotic value of about -0.4 cm/yr in 

the mantle below the asthenosphere. This is true both for 
/ 

slow and fast moving plates. For 150 My the asymptotic 
value of V drops below 0.1 cm/yr. 
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The upward- motion supplies the mass required for the 
accretion of the thickening plate. ’ One may rationalize 
the fact that the asymptotic value of v is essentially 
dependent upon age, but independent of the plate velocity, 
by writing, from the continuity equation, 

y a y 

/ /' ( 28 ) , 
o ^ o 


or, because x = (u • age) and v = 0, 

o o 


<3. f ,u ^ , 

•" d(age) u^^ ^ 


o o 


For a given age, the u/u^ profiles as plotted in Figures 2 
or 3 'are practically identical for slow and fast plates - 
Thus the above inregral is nearly independent of plate 
velocity and so is v^. This result follows immediately if 
one assumes that the plate thickness is a function of age 
only . Tie nexr. section will show that there are small 
departures from this approximation but that it is correct 
to within 20%.. The spread in values of v^ for a given 
age is also about 20%. 

A posteriori we may notice that the assiamption v « u 


is valid within the boundary layer for ages larger than 


10 My^ but would not hold much nearer to the ridge axis. 

We recall also that shear heating due to ^ — was neglected 
in (18) ; a posteriori we see this term to be small compared 


to 


8u 

9y* 
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Results for a Wet Olivine Mantle 

Calculations based on the olivine rheology given in 
the second column of Table 1 have been carried out to 
determine the effects of the presence of water in the mantle 
on the thermal and mechanical structure of the oceanic 
lithosphere and asthenosphere. Figure 7 compares the solu- 
tions for a wet and dry olivine mantle, both with V* = 

30 cm^/mole, T = 1200 and u — 10 cm/yr. The simple 
boundary layer cooling profile (conduction) for an age of 
150 My is shown for reference. The T and u profiles are 
only slightly different. The largest effect is at 150 My 
where the softer "wet" rheology reduces the amplitude of 
the viscous dissipation bump on the geotherm. 

The smaller value of E* in the wet olivine 

case yields values of the minimum viscosity ^ in the 
asthenosphere and. values of the shear stress x which are 
lower than those of corresponding dry olivine cases (see 
Table 4) . Again we see that the drop in viscosity- and 
shear stress is more marked at young age since at older 
ages shear heating can have some compensating feedback 
effect- on the -temperature. This feedback phenomenon was 
dominant in. our study of the continental asthenosphere 
(Froideva^uc and Schubert, 1975) . 

Solutxons for a Newtonian Rheology 

Six solutions for a Newtonian rheology with parameter 
values given in Table 1 are shown in Figures 8 (u and T) and 
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9 (ri) for a fast plate (u^ = 10 cm/yr) and a cold lower 
mantle (T = 1200 ^C) . The solid curves are for the New- 

CO 

3 

tonian I case with V* equal to either 11 or 30 cm /mole. 

The same convention as in previous figures identifies a 

particular value of V* . Recall that the shear zone is 

shallower for larger V*. The shaded regions again denote 

the range of solutions for intermediate values of V’', in 

the Newtonian I case. The dashed curves are the solutions 

for the Newtonian II set of parameter values. 

The temperature, velocity and viscosity profiXes can 

first be compared with the results for a dry olivine mantle 

shown at the bottoms of Figures 3 and 6, corresponding to 

the same values of u^ and T^. The viscous heating bumps 

-on the temperariire profiles are much more marked in the 

Newtonian cases and thus the depths of the shear zones 

are rather shallower. The relatively higher temperatures 

20 

lead to values of the viscosity minimixm between 10 and 
2X 

10 poises for rhe Newtonian I parameters. These values 
do not differ much from those of the corresponding dry 
non-Newtonian cases. Furthermore, the shear stress values 
are only a factor 1.5 to 2 larger than the corresponding 

T values for: the dry non-Newtonian cases at 10 My and 
become equivalent at 150 My (see Tables 3 and 5) . 

The comparison just discussed was between a rather 
readily deformable Newtonian rheology (Froidevaux and 
Schubert, 1975) and the relatively undeformable, dry 
olivine nonlinear case. We have seen that in the 
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presence of water the non-Newtonian rheology is much 
"softer". Therefore,' it may he more appropriate to com- 
pare the Newtonian 1 rheology with the wet olivine case, 
since both correspond to about the same value of the acti- 
vation energy E*. This comparison shows the non-Newtonian 
mechanism to be more efficient than the Newtonian I case, 
at least for the situation T = 1200 ?C and u =10 cm/yr. 

TO o 

The Newtonian II rheology leads to solutions with 
viscosities and shear stresses which are larger than those 
for the non-Newtonian cases, as can be seen in Figure 9 
and Table 5. The Newtonian II set of rheological parameters 
therefore also describes a mechanism of deformation v;hich 
is not. competitive with the nonlinear olivine rheology. 

This, remark applies again to a cold mantle (T^ = 1200 °C) 
and. .a fast plate. 

The Dominant Mechanism of Deformation 

The Newtonian law of deformation for solids is based 
on atomic diffusion across grains under the influence of 
applied stress (Nabarro, 1948) . In laboratory deformed 
olivine, this behavior has never been observed because the 
non-Newtonian process, based on the motion of dislocations 
(Weertman, 1970) , is much more efficient. Naturally de- 
formed peridot ite nodules exhibit sub-grain structures 
and dislocation densities similar to those of samples de- 
formed in the laboratory (Nicolas et al., 1973). Thus, 
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one might expect the non-Newtonian behavior to be dominant 
in the asthenosphere , wherein the peridot ite nodules are 
believed to originate. 

Our calculations can provide some insight to the 
question of what is the dominant mode of deformation in 
the mantle since we can compare values of shear stress 
for different sets of rheological parameters , other condi- 
tions being similar. However, there can be no definitive 
answer to this question since rheological parameter values 
for Newtonian deformation are purely speculative and we 
can only comment on the particular cases, Nex'?toni'an I and 
II , studied here . 

Figure 10 shows the values of shear stress at 1000 km 
from the ridge crest for a fast plate {u^ = 10 cm/yr) and 
for four different sets of rheological parameters as func- 
tions of the postulated temperature increment between the 
surface and the mantle below the asthenosphere - T^) . 

The shear stress for wet olivine {curve b) is always smaller 
than T for a dry olivine rheology (curve a) . Similarly, 
the Newtonian I (curve d) shear stress is always smaller 
than that of tdie Newtonian II case (curve c) . 

The important comparison one would like to make, how- 
ever , is between Newtonian and non-Newtonian cases for a 
dry or for a wet mantle. In the latter case, one compares 
curves d and b (if, in fact, the Newtonian I case represents 
the linear deformation of wet mantle rocks) . For 
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(T „ ip ) less than 1350 °C the non-Newtonian rheology 
yields smaller values of t than the Newtonian I rheology; 
for (T^ - T^) > 1350 the situation is reversed. Thus, 
at 1000 km from the ridge, i.e. for an age of 10 My, the 
dislocation mechanism is dominant for a relatively cold 
mantle, whereas for a relatively warm mantle the diffusion 
mechanism is more efficient. One can thus conceive a 
situation in which the deformation is essentially Newtonian 
near the ridge, where we shall see that stress values are 
small, and non-Newtonian under old oceans, where our cal- 
culations show stress values to be larger. 

We speculate that the Newtonian II rheology may be 
appropriate to a dry mantle. It may also be the case that 
the linear deformation of dry mantle rocks is even more 
difficult than that of the Newtonian II case. Figure 10 
shows that the dry non-Newtonian rheology yields values of 
T (curve a) smaller than those for the Newtonian II case 
(curve c) for (T - T ) < 1600 °C. Thus at 1000 km from 

CO O 

the ridge and for u^ = 10 cm/yr, the dislocation mechanism 
would be more efficient if the mantle is dry. 

The comparisons discussed above are illustrative of 
those that must be made at all distances from the ridge as 
a function of u^ and (T - T ) . Different mechanisms of 
deformation may be dominant in different locations rela- 
tive to ridge and trench and for different plate velocities 
and sub-asthenospheric temperatures . Non-Newtonian behavior 
certainly prevails in regions of large stress, whereas 
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Newtonian deformation takes over in regions of low stress 
On Figure 10 the cross-over occurs around 1 bar. For 
high enough temperatures at shallow depth;. (19) can stxll 
lead to xmportant deforraatxons even for low stress values 
The limitation on the present state of knowledge of the 
rheologxcal behavior of the mantle comes partly from the 
fact that laboratory experiments cannot be performed at 
low enough stress to measure the Newtonian parameters. 



40 


AGE DEPENDENCES OF THEBMAL AND MECHANICAL 
PROPERTIES OF THE LITHOSPHERE AND ASTHENOSPHERE 


Lithosphere Thickness 

In discussing the thickness of the lithosphere and 
how this thickness increases with age, one must first be 
precise about the definition of this quantity. Since the 
term lithosphere is generally understood to refer to the 
essentially rigid plates at the Earth’s surface, a mechan- 
ical definition of lithosphere thickness seems appropriate. 

We define lithosphere thickness to be the depth at which 
u/u^ = 0.95; above this depth there is essentially no 
de form ation , while below it there is a region of strong 
-shear, the asthenosphere. Others, because of either the 

geophysical data they consider (seismic or heat flow) or 

* 

limitations of the theoretical model they use to describe 
the Lithosphere (the mechanical structure of the lithosphere 
has not heretofore been modelled) have employed the term 
lithospheric thickness to indicate the depth of a particular 
isotherm or a characteristic depth of a seismic wave velocity 
model. 

Figure 11 shows the thickness of the rigid lithosphere, 

h 

defined above, as a function of (age) and age, for differ- 
ent rheologies , rheologxcal parameter values , plate velocities 
and for a deep temperature T^ = 1200 °C. The solid curves 
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are for the dry non-Newtonian rheology; for the upper two 

3 

solid curves V* = 11 cm /mole and for the lower two solid 

3 

curves V* = 30 cm /mole. Plate velocity distinguishes 
between the two solid curves for each value of V*; the 
one with the larger thickness is for u^ = 2 cra/yr, the 
other is for 10 cm/yr. The dashed curves are- for 

the Newtonian I rheology and u^ = 10 cm/yr; the upper 

3 

dashed curve is for V* = 11 cm /mole, the lower one is 

for V* = 30 cm^/mole. Finally, the long- short dashed 

curve is for the Newtonian II rheology and u^ = 10 cm/yr. 

The lithosphere is observed to thicken with age, but 

r . 

it grows more slowly than (age)'^. As we noted earlier, the 
lithosphere is thicker the smaller the activation volume 
and the smaller nhe plate velocity, other quantities being 
the same. In general, the rate of growth of the lithosphere 
decreases with age, an effect which appears to be more pro- 
nounced under conditions wherein viscous dissipation con- 
tributes significantly to the therrao-mechanical structure 
(thus, for 1600 the effect would be less apparent) . 

For example in the case of the dry non-Newtonian rheology 
with V* = 30 cm?/mole, T^ = 1200 and u^ = 10 cm/yr, the 
lithosphere thickness at an age of 10 My is 40 km. If 

j- 

this lithosphere thickened with (age) ^ it would be about 
155 km -chick at 150 My; however, according to Figure 11 it 
is only 100 km thick at 150 My. The thickening of the 
lithosphere with age is mainly the result of boundary layer 
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cooling. Viscous heating, especially at older ages tends 
to reduce the rate of growth of the boundary layer or 
lithosphere . 

Although, as discussed above, we cannot directly 
compare our determinations of lithosphere thickness based 
on the mechanical definition with that inferred from model 
inversions of surface wave phase velocities (Leeds et al.,1974; 
Forsyth, - 1975 t Leeds, 1975), a qualitative comparison may be worth 
noting. Leeds (1975) concludes that the lithosphere con- 
tinues to thicken even 100 My after formation, in agreement 
with the results of our study. A quantitative comparison 
of the rate of thickening is unwarranted in view of the 
differences in definition of Ixthosphere thickness and 
-the dependence of the seismically inferred thickness on 
the model used in the inversions. Another piece of 
evidence for a continuous thickening of the oceanic 
lithosphere may be rhe fact that old plates require a 
longer time to be thermally equilibrated after they plunge 
into the mantle at the trenches (Deffeyes, 1972) . 

Our results are to be contrasted with those of 
McKenzie (19ff7) and Parsons and Sclater (1976) who assume, 
a priori , a lithosphere of constant thickness and 
Oldenburg (1975) who finds that the lithosphere thickens 
proportionally with (age) 
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Heat Flow 

Figures 12 and 13 show the surface heat flux in 
2 

HFU (1 ucal/cm -sec) as a function of age in My. Also 

included in the figures are curves of shear stress as a 

function of age; these will be discussed in a separate 

sub-section. The heat flow curves show values of a 

* o 

which decrease with age. Figure 12 presents the heat flow 
for the non-Newtonian rheology; cases a and b are dry 

3 

olivine and wet olivine, respectively, with v* = 30 cm /mole. 
The curves in the upper part of the figure are for u^ = 

2 cm/yr and those in the lower part have u^ = 10 cm/yr. 

Heat flux results shown in Figure 13 are ;based on the 

Newronian rheologies, the solid curves for the Newtonian I 

/ 

.param^er values and the dashed curve for the Newtonian II 
case^ The numbers 30 and 11 labelliiig the solid curves 
in Figure 13 refer ’to -the values of, activation volume in 
cm /mole. The_ deep temperature for the results of both 
Figures 12 and 13 was again 1200 °C. Finally, the dotted 
lines in the figures show the heat flow as a function 
of age for the simple boundary layer cooling situation, 
wherein the surface heat flux decreases with the reciprocal 
of (age) The dotted lines therefore have slopes of ~h. 

The solutions of Figures 12 and 13, which include 
theormo-mechanical coupling and viscous dissipation, show 
nouiceable departures from the inverse square root 
dependence on age. The departure of q^ from- the inverse 
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square root dependence on age becomes more pronounced 
with increasing age; all the curves have an upward curva- 
ture and show a tendency to flatten at old ages. There 
is a dxfference between the of the simple boundary layer 
coolxng solution and the solutions of Figures 12 and 13 even 
at young ages. Particularly apparent is the excess heat 
flux at 10 My and older associated with viscous dissipa- 
tion in two of the Newtonian cases of Figure 13 and in the 
non-Newtonian case of Figure 12, ~ 10 cm/yr {curves a 

and b in rhe lower part of Fxgure 12) . For the slowly 
moving plate (u^ = 2 cm/yr, upper part of Figure 12) , shear 
heating effects for the wet olivine rheology are noticeable 
only axter about 100 My, while the dry olivine rheology 
leads to appreciable viscous dissipation beyond about 20 My. 

The heat generated by viscous dissipation adds to 
than, of boundary layer cooling and tends to displace the 
curve upward from rhe straight line associated with 
houndary layer cooling, and because shear heating becomes 
more important with age, it tends to flatten the q^ curve 
with increasing age. A comparison of the Newtonian curves in 
Fx-gure 13 shows that the larger V*, the larger the contribution 
of vxscous dissipation to q^ and the greater the flattening 
of the q^ vs. age curve at older ages, other quantities 
being equal. Fxgure 12 shows that there is a modest increase 
in q^ for hxgher values of plate velocity. 
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Sclafcer et al. (1976) have reconsidered the age dependence 

of heat flow averages for ohe North Pacific. Their heat flow 

averages included data only from sites with a thick layer of 

"impermeable” sediment, since it is generally believed that hydro- 

thermal circulation can transport a significant: amount of heat 

(Elder, 1965) whereas the measurements only reflect that portion 

of the total heat being conducted through the sediments (Lister, 

1972, 1974). Figure 14 compares these average heat flows, shown 

by the shaded rectangles, with three theoretical surface heat 

flux curves derived from our model. The vertical and horizontal 

dimensions of each box are estimates of the uncertainties in 

heat flow and age, respectively (Sclater et al., 1976). The 

filled circles wirh associated error bars are data for the 

-Nornh Pacific from Sclater and Francheteau (1970) . The solid 

CU3TVS is for a wer olivine rheology with an activation volume 

of 11 cm mole “ and = 1600 °C; the long and short dashed 

3' -1 3 —1 

curves are foxr V* — 30 cm mole and 11 cm mole , respectively, 
with a dry olivine rheology and T^ = 1200 °C. The plate velocity 
in all three cases is 10 cmyr appropriate for the fast Pacific 
plate - 

The averages of Sclater et al. (1976) are consistent with 

— 

a heat flow proportional to (age) ^ for most of the ocean floor 
and suggest a flattening of the heat flow in the ocean floor 
older than about 100 My. These heat flow averages are in quali- 
tative agreement with the results of our model. With T = 1200 

OO * 

viscous heating produces a flattening of the theoretical heat 
flow curves for the older ocean floor. For T = 1600 ^C, shear 



46 


heating does not contribute signifxcantly to the theoretical 

surface heat flux. One might argue that available heat flow 

values for old oceans could be used to determine what T or 

00 

what rheological parameters are appropriate for the Earth. 
Figure 14 shows, however, that oceanic heat flow measurements 
are not yet accurate enough to discriminate among the various 
cases. In particular, we do not feel that can presently 
be inferred by this method. Additional factors not considered 
in our present model, such as mantle radioactivity, further 
complicate the task of choosing preferred values of or 
rheological parameters. 

Shear Stress 

_ The dependence of shear stress x on age of the ocean 
floor is also shown in Figures 12 and 13; parameter values 
for the different curves follow the same convention as given 
previously for the heat flux values. Shear stress is seen, 
in general, to increase with age of the ocean floor; for the 
Newtonian IX rheology, however, there is a slight decrease 
of T with age (Figure 13) . A larger activation volume implies 
a larger x, other parameters being equal. The shear stress 
for dry olivine is larger than that of the wet olivine for 
the same value of V*. 

The magnitude of x lies between about 10 and 100 bars for 
the cases shown in Figures 12 and 13 (T^ = 1200 °C). For very 
young ages, higher T^ values and for more readily defoarraable 
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rocks, e.g. a wet olxvine with V* = 11 cut /mole, x can be as 

small as a fraction of a bar to several ' bars . We have -also 

seen in Figure 10 that higher sub-asthenospheric temperatures 

lead to lower values of shear stress. Solomon et al. (1975) 

have inferred shear stresse-s of a few bars at the base of 

oceanic plates from their studies of intraplate stresses. 

One may integrate the shear stress vs.' age relationships 

given in the figures to obtain the total drag force per unit 

length of ridge acting on a lithospheric plate. For both the 

Pacific and Atlantic plates the total drag force per unit length 

X6 17‘ 

of ridge would be of the order of 10 -10 dynes/cm if condi- 

tions, are such that the x values of Figures 12 and 13 apply. If' 

\ 

this had, to be balanced by a pull from the downgoing slab at the 

/ 

tresich.,. it would imply a downward force per unit depth and per unit 
length, of trench of 1-10 kbars . This may be somewhat too high al- 
though,' Schubert et al. (1975) have found that the downward force on 

i 

the slab due to thermal contraction land phase boundary elevation may 
be several kbars per unit length of trench and per unit depuh. The 
drag force may be one to two orders of magnitude smaller, especially 
for higher values of T^. Froidevaux and Schubert ’(1975) found shear 
stresses beneath continents to be of the order of 200 bars for tem- 
perarures at 400 km comparable with T^ values considered here. The 
addinion of radioactive heat sources in the continental mantle re- 
duced the shear stress to 40 bars.. It is thus fair to say, by com- 
parison, that the stresses under oceanic plates are smaller by one 
order of magnitude or more. It is quite possible however that 
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different rheological properties prevaxl under the two 
types of provxnces. Our conclusion is in agreement with 
that of Forsyth and Uyeda (1975) based on their study of 
the relative forces acting on plates in dynamic equilibrium. 
They found that the drag resistxng plate motion 

is stronger under the contxnents than under the oceans. 

The type of data required to uhrow some light on this 
problem is a detailed comparative study of the disloca- 
tion structure and density in the sheared peridotite 
nodules from oceanic and continental origin. 

We have noted earlier that dr/dx must be small com- 
pared with pg if the mathematical formiulation of our 
model, is valid. Prom Fxgures 12 and 13 one finds that 

t 

- dx/dx- is at most about 0.5 dynes/cm and is usually sig- 

3 3 

nj.rxcantly less, whereas pg is about 3x10 dynes/cm . 

Thus our approximation dx/dx << pg is an excellent one; 

i 

it says that the pressure is purely lithostatic. 

The net state of stress in the lithosphere is a 
complex superposition of thermal and mechanical stresses 
obviously not accounted for by our model. This precludes 
a direct comparison of intraplate stresses inferred from 
earthquake source mechanisms with the t values computed 
here. Since intraplate stresses inferred from earthquake 
source mechanisms and xn-sxtu measurements reflect the 
resultant, in the elastic part of the lithosphere, of 
various large scale forces acting at plate boundaries 
(Sbar and Sykes, 1973) and of more local stresses, they 
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cannot directly compare with the shear stress values of 
our model. Our shear stress values . represent the viscous 
drag acting beneath the plates . 

Bathymetry 

As discussed in the section on the mathematical 
formulation of our model, the themal contraction and densi- 
fication of the lithosphere as it moves away from the ridge 
and cools, can be interpreted as giving rise to the topography 
of the ocean floor via the mechanism of isostatic compensa- 
tion. Using equation (4) , we have calculated the depth of 
the ocean floor as a function of age for the dry non- 
Newtonian case with = 1200 °C, u^ = 10 cui/yr and V* = 

3 

11 and 30 cm /mole. The results are shown in Figure 15 

where the depth of the ocean is plotted as a function of 

(age)^; the upper part of the figure is for an expansion 

-5 -1 

coefficient a = 3x10 K and the lower part is for a ~ 

-5 -1 

4x10 K . The shaded region again represents the range 
of topographic solutions for intermediate values of V*. 

The figure also includes the topography associated with 
the simple boundary layer cooling solution computed on 
the basis o:f (16) (the dashed curves) . Finally, bathyme- 
tric data for the Pacific ocean, taken from the summary of 
Parsons and Sclarer (1976) , are shown on the figure (open 
circles are based on depths from topographic profiles, 
closed circles are based on depths from Joides sites) . 

Since depth increases as (age)"^ in the simple boundary 
cooling, solution (equation (16)), the slopes of the straight 
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dashed lines indicate this functional dependence. Note that 
the slope depends on a and thus is different in the upper 

I 

and lower parts of the figure. It is clear that the data 
depart from the straight line for ages greater than 80-100 My. 

The ocean floor beyond 100 My is not as deep as would be pre- 
dicted by simple boundary layer cooling. It is also' obvious 
in the figure that our thermo-mechanical model of the oceanic 
lithosphere and asthenosphere yields depths which depart sig- 
nificantly from the (age) law^ showing considerable flattening 

as a function of age . In the upper part of the figure, for 
-5 -1 

a = 3x10 K / it is seen that no reasonable value of V* would 

yield topography in good agreement with the data while in the 

“5 “X 

lower part of the figure, for a = 4x10 K , one can see that 

3 

a Y* somewhat larger than 11 cm /mole would adequately match 
the observations. These statements are true for the boundary 
condition T^ = 1200 °C. For higher values of T^ shear heating 
effects are less pronounced and a higher V* value and/or a smaller 
a value are more appropriate to fit the data. Such trade-offs of 
parameters are often characteristic of geophysical modelling. 

The flattening of ocean floor topography with age is 
a consequence of the viscous dissipation in our model. 

Frictional heatring in the shear zone that constitutes the 
asthenosphere offsets the conductive cooling of the lith- 
osphere as it moves av7ay from the ridge. The shear heat- 
ing tends to maintain the material at a higher temperature 
as a function of age than if the material were simply 
cooling by conduction. Thus less thermal contraction 
occurs and the lithosphere remains lighter with age; 
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the depth of the ocean floor therefore need not .increase 
as rapidly as (age) ^ to maintain a state of isostatic com- 
pensation. 

The model of Parsons and Sclater (1976) also shows an 
ocean floor topography that tends to flatten with age'. 
However, their result is a consequence of an a-priori speci- 
fied feature of their model, namely that the lithosphere 
is of constant thickness for all ages. There is no physi- 
cal mechanism in their model, such as the shear heating 
in ours, which is responsible for the flattening of the 
topography. We recall that our results, discussed at the 
beginning of this section, show that the lithosphere con- 
tinues uo increase in thickness even at old ages. 

Seismic Velocities 

Our model is that of a pure olivine mantle. This 

i 

f 

approximation is justified by the fact that olivine is 
the- main constituent of the upper mantle, and that obser- 
vaxions of peridotite nodules (Nicolas et al. , 1973} show 
thar olivine is more abundant and deforms more readily than 
the pyroxenes and other minerals contained therein. How- 
ever, an investigation of the elastic properties of the 
oceanic mantle certainly requires a more adequate petrologi- 
cal model and the inclusion of possible phase changes. In 
particular, it should be decided whether or not the asthen- 
osphere is a zone of partial melting. Our model does not 
require partial melting for a zone of mechanical decoupl- 
ing, since solid state deformation yields sufficiently low 
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vxscosities. On the other hand, it is known from the 
study of metallxc phases (Auten and Gordon, 1975) that a 
small amount of partial melting does not substantially 
alter the subsolidus deformation law, since the grains 
must still be deformed in their bulk to allow shear on a 
large scale. This is also true for rocks (Arzi, 1972) . 

Hence, we have not, till now, been concerned with whether 
the computed geo therms cross the solidus of mantle material. 

Figure 16, representing seismic velocities vs. depth 
for ages of 10 My and 150 My, has been drawn on the basis 
of elastic parameters for pure dry olxvine (Kumazawa and 
Anderson, 1969) . The solid lines with the most pronounced 
velocity minima are for = 1600 the other solid 

linee are for T = 1200 and the area in between is 

' oo 

shaded; both cases are for a slow plate (u^ = 2 cm/yr) . 

The dashed solutions are for u^ = 10 cm/yr, = 1200 
3 1 

and = 30 cm /mole, i.e. for a geotherm exhibiting a well 
developed shear heating bump. At 10 My this case falls on 
the solid line for = 1200 °C. For their prefered con- 
tinental geotherm under old shields, Proidevaux and Schubert 
(X975) obtained a flat seismic velocity profile in the 
asthenosphere. Here, on the contrary, a marked minimum 
develops around 50 km depth for an age of 10 My and below 
100 km depth for an age of 150 My. This contrast between 
oceanic and continental upper mantle seismic velocity pro- 
files is in good agreement with seismologxcal evidence 
(Knopoff, 1972; Press, 1972; Jordan, 1975). Fur thermore , the amplitudi 
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of this seismic velocity drop, as seen from Figure 16, is 
of the right magnitude, between 0-2* and 0.45 km/sec, for 
shear waves. The presence of a molten phase would of 
course amplify these drops", perhaps by an embarrassingly 
large amount. 

Studies of Rayleigh waves (Leeds, 1975) are usually 
analyzed by assuming a constant value of v in the top part 
of the upper mantle, the seismologist’s lithosphere, and a 
value 0.5 km/sec lower in the low velocity channel, the 
seismologist’s asthenosphere. Our model suggests that the 
velocity drop, obtained on the basis of pressure and tem- 
perature effects alone for an unsophisticated petrology, is 
smaller for older ages, as can be seen in Figure 16. This 
-conclusion would still be qualitatively correct if partial 
melting occurred in the shear zone. The flattening of 
the geotherm with age implies that partial melting would 
be more pronoimced, or restricted, to young ages. For 
example,^ Figure 17 shows that the temperature profiles 
of the purely thermal solutions of (7) for T^ = 1200 °C 
do not 'intersect the solidus curve for a dry mantle 

(curve a) . in the presence of 
hydrared mmerals (curve b) this conclusion still holds. 

For Jiigher vaXues of T^ the temperature of the solidus is 

^ i « 

exceeded, in particular for young ages. Curve 'c', in ‘Figure 17 
is the solxdus for a completely wet mantle. The solid lines 
a, b and c are from Kushiro et al. (1968) . The dashed 
solidus is that given by Ito and Kennedy (1967) for dry 
peridotite- 
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These remarks are not applicable to the region of 
the ridge crest, since our model is for ages larger than 
10 My. At the ridge crest, the high attenuation of 
seismic waves testifies to the presence of a fraction of 
melt {Solomon, 1973) , also required to generate the 
tholeitic basalts forming the new oceanic crust. What 
is strongly needed is a series of seismological experi- 
ments parallel to the ridge axis with ocean-bottom seis- 
mometers. This should help to settle the problem of the 
extent of partial melting in the oceanic asthenosphere. 

The low velocity channel corresponding to strong 
temperature gradients in the oceanic lithosphere does 
not coincide, in depth, with the shear zone. This can 
-be seen by comparing the depths of the seismic velocity 
mirrina. in Figure 16 with the depths of the viscosity minima 
in T-lgures 3 or 6. The seismologist's low velocity zone 
IS Trh.us probably shallower than the*'- asthenosphere defined 
as the zone of decoupling of the plate from the mantle 
below. 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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CONCLUDING REMARKS 


Previous models of the oceanic lithosphere have all 
been purely thermal, x.e. the velocity fields are not 
derivable therefrom in a self-consistent manner owxng to the 
coupling of thermal and mechanical fields via the depend- 
ence of viscosity on temperature and via viscous dissipation. 
The main emphasis in this paper has been on the development 
of a coupled thermoraechanical subsolidus model of the oceanic 
lithosphere and asthenosphere and on the exploration of the 
consequences thereof. "For example, what departures from 
the simple (age) “ dependences for such observables as ocean 
floor depth result from thermomechanical coupling? Also, 
-what are the velocity profiles beneath the oceans and what 
are the' mechanical lithosphere and asthenosphere thicknesses 
derived therexxom? 

We have found that simple boundary layer cooling with 
its direct and inverse proportionalities of ocean bottom 
depth and "ocean floor heat flux on (age) is predominant at 
young ages (not too near the ridge however) . At older ages, 
effects of viscous dissipation become more important, since 
lithospheric cooling and thickening tends to place the 
asthenosphere at: greater depths where the pressure-dependence 
of the viscosity leads to higher values of viscosity in the 
asthenosphere. Shear heating causes the ocean to deepen 
less rapidly than (age) the lithosphere to thicken more 

p 

slowly than (age) and the surface heat flow to decrease 
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less rapidly than (age) The magnitudes,© f the shear 

heating effects depend not only on age, but also on the 
plate velocity and the activation volume of the deforma- 
tion process. 

Shear stresses and values of the minimum viscosity 
in the asthenosphere generally increase as the age of the 
ocean floor and are substantially smaller under the oceans 
than under the continents , again mainly because of the 
greater depth of the continental asthenosphere coupled 
with the pressure-dependence of the viscosity. It is im- 
portant to remember that the depth of the asthenosphere is 
not imposed a priori but is calculated. Thus the drag 
forces beneath plates with continents are larger than 
those. beneath oceanic plates. 

Because of our emphasis on studying the physical 
effects of thermomechanical coupling on lithosphere and 
asthenosphere structure, we have not attempted to adjust 
unknown parameters such as V* and T^ to obtain the best 
quantitative fit of our model results with observations of 
such quantities as heat flux and topography. Nor 

have we tried to distinguish which of the rheological laws 
provide the best fit to such data. Our qualitative com- 
parisons with data leave no doubt that parameters can be 
adjusted to have physically reasonable values which pro- 
vide good quantitative fits to heat flow and topographic 
data, for example. In fact, with V*, T^ and the rheological 
law to vary, in addition to other parameters (e.g. a in 



57 


the case of ocean floor topography) , it is probably true that 
many different combinations of parameter values will pro- 
vide adequate fits to data. 

More importantly, a number of physical effects not in- 
cluded in our present study should be investigated before 
attempting to quantitatively fit observations. For example, 
the influences of radioactive heat sources in the oceanic 
upper mantle and crust and a nonzero heat flow from great 
depth need to be determined. We have already seen in our 
study of the continental asthenosphere that the addition 
of radioactivity could significantly reduce the shear 
stress. 

We need also to investigate the effdcts of a possible 
shaircw return flow in the asthenosphere forced by a non- 
zero adverse pressure gradient. The approach of the pre-. 

I 4 

sent' paper can readily be modified ,to include this effect 

\ 

thereby extending the work of Schubert and Turqotte (1972) . 

In particular', it is important to d^emineif a shallow return flow 
rs compatible with gravimetric and bathymetric data on the 
basis of the more complete approach developed here. As 
regards shear heating effects and associated departures 
from (age)-^ laws it is clear that they are more important 
in the return flow case. 

Our present study of the oceanic lithosphere and 
asthenosphere, like our previous investigation of the 


omginalwgbb 
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continental asthenosphere, emphasizes the importance to 
geodynamic modelling of acquiring more laboratory data 
on the thermal and mechanical propertxes of mxnerals pre- 
sent in the mantle. Particularly important in this connec- 
tion are determinations of the activation volume for olivine 
In terms of field measurements/ seismic observations which 
yxeld velocity profiles and attenuation data as a function 
of ocean floor age are vital for further progress in our 
understanding of the structure of the oceanic lithosphere 
and asthenosphere . It is equally hoped that studxes of 
phase equilibria and deformation patterns in peridotite 
nodules found in kimberlites (Africa/ Siberia) and in 

basal-tic lava flows (volcanic islands and. continents) will 

< 

soon provide more constraints on the temperature and vel- 
ocity profiles , as well as on the magnitudes of shear stress 
A major question which such observations could help answer 

I 

is the extent to which the asthox);5)here may be partially 
melted away from the ridge. 



59 


ACKNOWLEDGMENT 


This work was partly supported by the Earrh Sciences 
Section, National Science Foundatxon, NSF Grant GA 40749; by 
the National Aeronautics and Space Administration, NSG 7002; 
and by the University of Paris-Sud, Orsay where Gerald 
Schubert was a Visiting Pxofessor. 



60 


REFERENCES 

Anderson, D. L* and C. Sammis, Partial melting in the upper 
mantle, Phys. of Earth and Planet. Interiors , 41-50, 
1970. 

Arzi, A. A. , Experimental study of partial melting in natural 
rocks and subsequent creep under low stress, Trans. Am. 
Geophy s . Union , 53 , 513, 1972. 

Auten, T. A. and R, B. Gordon, Compressive creep rate of par- 
tially melted AlGa alloys, Metall . Trans . , 6A , 584-586. 
1975. 

Clark, S. P., Jr., editor. Handbook of Physical Constants , 
Geological Society of America, Memoir 97, p. 92, 1966. 

navis, E. E. and C. R. B. Lister, Fundamentals of ridge 

crest topography, Earth Planet. Sci. Lett. ,'' 21, 405- 
413,. 1974. 

Deffeyes, K. S., PXxime convection with upper mantle temper- 
ature inversion. Nature , 240 , 539-544, 1972. 

Elder, J. W. , Physical processes in geothermal areas, in 

Terrestrial Heat Flow , edited by W. H. K. Lee, pp. 211- 
239, American Geophysical Union, Washington, D. C., 1965. 

Elsasser:, W. M. , Convection and stress propagation in the 
upper" mantle., in The Application of Modern Physics 
to the Earth and Planetary Interiors , edited by S. K. 
Runcorn, pp. 223-246, Wiley-Interscience , New York, 


1969. 



61 


Forsyth, D- , The early structural evolution and anisotropy 
of the oceanic upper mantle, Geophys. J, Roy. Astron. 
Soc. , 103-162, 1975. 

Forsyth, D. and S. Uyeda, On the relative importance of the 
driving forces of plate motion, Geophys. J. Roy. 

Astron. Soc.' , 43 , 163-200, 1975. 

Froidevaux, C. and G, Schubert, Plate motion and structure 

of the continental asthenosphere: A realistic model of' 

the upper mantle, J. Geophys. Res., 80-, 2553-2564, 1975 

Gueguen, Y. and J. M. Mercier, High attenuation and’ low 

velocity zone, Phys. Earth and Planet. Interiors , l_f 
2L39, 1972. 

I to, K.. and G. C. Kennedy, Melting and phase relations in 

a natural peridotite to 40 kilobars. Am. J. Sci . , 265 , 
5X9-538, 1967. 

Jordan, T. H. , The continental tectosphere. Rev. Geophys. 

I - — ' 

Space Phys. , 13, 1-12, 1975. ' 

\ 

KauseX, E. G. , A. R. Leeds and L. Knopoff, Variations of 

Jlayleigh wave phase velocities across the Pacific ocean. 
Science , 186 , 139-141, 1974. 

Knopoff, L. , Observation and inversion of surface v?ave disper 
sion, Tectonophysics , 13 , 497-519, 1972. 

Kohls tedt, D. L=. and C. Goetze, Low-stress and high tem- 
perature creep in olivine single crystals, J. Geophys. 
Res. , 7^, 2045-2051, 1974. 

Kohlstedt, D. L. , C. Goetze and W. B. Durham, Experimental 
deformation of single crystal olivine with applica- 



62 


tion to flow in the nyantle, in Petrophysics : The 

Physics and Chejnistry of Minerals and Rocks , ed . 

S. K. Runcorn, John Wiley & Sons, Ltd., London, 1976. 

Kumazawa, N. and 0. L. Anderson, Elastic moduli, pressure 
derivatives and temperature derivatives of single 
crystal olivine and single crystal forsterite, J. 

Geophys. Res. , 75 , 5961-5972, 1969. 

Kushiro, I., Y. Syono and S. Akimoto, Melting of a peridotite 
nodule at high pressures and high water pressures, 

J. Geophys. Res. , 73 , 6023-6029, 1968. 

Langseth, M, G. , Jr., X. LePichon and M. Ewing, Crustal 

structure of mid-ocean. ridges, 5. Heat flow through 
the Atlantic ocean floor and convection currents, 

CT. Geophys. Res. , 71 , 5321-5355, 1966. 

Leeds:, A. R. , Lithospheric thickness in the western Pacific, 
Rhys. Earth Planet. Interior , 11 , 61-64, 1975. 

Leeds., A. R. , L. Knopoff and E. G. Kausel, Variations of 
upper mantle structure under the Pacific ocean. 

Science , 186 , 141-143, 1974. 

Lisuer, C. R- B. , On the thermal balance of a mid^ocean ridge, 
Geophys. J. Roy. Astr. Soc. , 26 , 515-535, 1972. 

Lister, C. R. B. , On the penetration of water into hot 

rock, Geophys. J. Roy. Astr. Soc. , 39 , 465-509, 1974. 

McKenzie, D. P., Some remarks on heat flow and gravity anom- 
alies, J. Geophys. Res. , 72 , 6261-6273, 1967. 

Nabarro, E. R. N. , Deformation of crystals by the motion of 
single ions, in Conference on the Strength of Solids, 



63 


pp, 75-90, Physical Society, London, 1948. 

Nicolas, H. , F. Boudier and A. M. Boullxer, Mechanisms of 

flow in naturally and experxmentally deformed perxdotites, 
2\m. J. Sci. , 273 , 853-876, 1973. 

Oldenburg, D. W. , A physical model for the creation of the 
lithosphere, Geophys. J. Roy. Astr. Soc. , 43 , 426-451, 
1975. 

Parker, R. L. and D. W. Oldenburg, Thermal model of ocean 
ridges, Nature , 242 , 137-139, 1973. 

Parsons, B. and J. G. Sclater, An analysis of the variation 
of ocean floor heat flow and bathymetry with age, s\ib- 
mitted to J. Geophys. Res. , 1976. 

Post, R. L. , The flow laws of Mount Burnett dunite, Ph.D. 
-thesis, Unlv. of Calif., Los Angeles, 1973. 

Press, F., The earth's interior as inferred from a famxly of 
of models, in The Nature of the Solid Earth , edited by 
E. C. Robefrson, pp. 147-171, McGraw-Hill, New York, 

1972. 

Sbar, M. L. and L. R. Sykes, Contemporary compressive stress 
and. seismicity in Eastern North America; an example of 
intraplate tectonxcs, Bull, Geol. Soc. Am. , 84 , 

18&1-1882, X973. 

Schatz, J. F.. and G. Simmons, Thermal conductivity of earth 
materials at high temperatures, »T. Geophys. Res. , 77 , 
696&-6983, 1972. 

Schubert, G. and D. L. Turcot te. One-dimensional model of 
shallow mantle convection, J. Geophys. Res., 77, 945- 



64 


951, 1972. 

Schubert, G. , D. A._ Yuen, and D. L. Turcotte^ Role of phase 

transitions in a dynamic mantle , Geophys . J. Roy. Astron. 
Soc. , 705-735, 1975. 

Sclater, J. G. and J. Francheteau, The implication of ter- 
restrial heat flow observations on current tectonic and 
geochemical models of the crust and upper mantle of the 
earth, Geophys. J, Roy. Astro. Soc., 2^, 509-537, 1970. 

Sclater, J. G., J. Crowe and R. N. Anderson, On the reliability 
of oceanic heat flow averages, J. Geophys. Res., in press, 
1976. 

Solomon, S. C. , Shear wave attenuation and melting beneath the 
mid-Atlantic ridge, g. Geophys. Res. , 78 , 6044-6059, 1973. 

Solomon, S. C. and N. H. Sleep, Some simple physical models for 

absolute plate motion, J. Geophys. Res. , 79 , 2557-2567, 1974 

Solomon, S. C- , N. H. Sleep and R. M. Richardson, On the forces 
driving plate tectonics; Inferences from absolute plate 
velocities and intraplate stress, Geophys. J. Roy. Astr. 

Soc . , in. pr-ess, 1975. 

Turcotte, D. L. and E. R. Oxburgh, Finite amplitude convective 

cells and concinental drift, J. Fluid Mech. , 28 , 29-42, 1967 

Turcotte-, D. X.. and E. R. Oxburgh, Mantle convection and the 
new global tectonics, Ann. Rev. Fluid Mech. , 4, 33-68, 

1972. 

Weertman, J. , The creep strength of the earth’s mantle. Rev, 
Geophys . Space Phys . , 8, 145-168, 1970. 



RY NON-NEWTONIAN 
(OLIVINE) 


kcal 

mole 


WET N0N“NEWT0NIAN 
(OLIVINEl) 


NEWTONIAN I 


kcal 

mole 


kcal 

mole 


NEWTONIAN II 



kcal 

mole 


« 11-30 


mole 


,=6.'45x10'13 b,= 5 . 41x10^^5 omSg eA g cnigeoK 

3 3, 3 3 1 gm 


= 11-30 


mole 


= 11-30 


mole 



E -l.sixio-^ 25S|2l< 
X 5^ 


ohlstedt and Goetze 
(1974) and 
Kohlstedt et al. 
(1975) 



Froidevaux and 


Schubert (1975) 


Schubert and 


Turcotte (1972) 


‘Table It Rheological Parameter Values 

















V* (■'^2^") 
mole 


o 

(cm/yr) 



11 

1200 

' 30 

1 

1 

11 i 

1600 

30 

°c‘ 

3.16* 

3.16* 

4.23* 

4.23* 

3.17 

" 

3.23 

1 

4.23* 

4.23* 

0.82* 

1 

1 0.93 

1.09* 

1.09* 

0.'8S ■ 

1 1.08 
1 

1.09* 

1.13 


1 

! 




Table 2. Surface heat flow values for a dry non-Newtonian 
rheology. Here is defined positive upward and 
is given in HFU (ycal/cm^* sec) , The asterisk denotes 
that the value is identical to that of simple boundary 
layer cooling. 
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Comparison of viscosity minimum nmin and shear 
stress T values for wet and dry olivine rheologies 
with T =1200°C, u =10 cm/yr and V*=30 cm^/mole. 
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46.2 
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91.4 


Table 5. Values of shear stres's (bars) for the Newtonian solutions 
(T =1200°C, u =10 cm/yr) . 
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FIGURE CAPTIONS 

Figure 1. Sketch of the model showing the thickening 

oceanic lithosphere and the asthenosphere be- 
neath. The left hand side of the diagram lists 
the quantities that are inputs to the model; 
the right hand side of the figure shows those 
quantities which are outputs from the calcula- 
tions. Note the two-dimensional flow in the 
asthenosphere; the vertical motion supplies the 
material to the growing lithosphere. T is temper- 
ature, u, V are horizontal and vertical veloci- 
ties, T is shear stress. There are no sharp 
boundaries to the asthenosphere. The boundaries 
shown here have artistic value only. 

Figure 2. Temperature T and horizontal velocity u pro- 
files for a -dry olivine rheology as defined in 
Table 1, column 1, T =^1600 °C and u = 2 and 
10 cm/yr for the top and bottom parts, respec- 
tively. Solutions are plotted for two extreme 

3 

values of V*, 11 and 30 cm /mole, and for two 
ages, 10 and 150 My. Temperature profiles are 
insensitive, at the accuracy of our graph, to 
the value of V*. This is not so for the velocity 
curves. For a given age and u^ the lower u vs. 

3 

depth curve always corresponds to V* == 11 cm /mole, 

3 

the upper one to V* = 30 cm /mole. The shaded 
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area between two such curves represents the 
range of solutions for intermediate values of 
V*. The T profiles obtained for the simple 
boundary layer cooling problem are identical 
with the two profiles seen at the top of this 
figure. For u^ = 10 cm/yr and an age of 150 My 
the above T vs. depth curve differs significantly 
from the boundary layer cooling solution. 

Figure 3. Temperature T and horizontal velocity u pro- 
files for a dry olivine rheology. The format 
is identical to that of Figure 2. Here = 

1200 °C. For an age of 150 My the temperature 

profiles are now split. The upper curve is 
3 

for V* = 30 cm /mole/ and the lower one is for 
V* = 11 cm^/mole. The splitting of the u pro- 
files is as in Figure 2 and the shading again 
indicates the range of solutions for intermedi- 
ate V* values. The dotted curves represent 
the temperature profiles for simple boundary 
layer cooling with variable thermal conductivity. 
In- the top part (slow plate) these solutions 
are- practically identical to the 10 My and 

3 

150 My, V* == 11 cm /mole, curves, respectively. 
Figure 4-. Surface heat flow vs . the temperature differ- 
ence between the deep mantle and the surface 
(T - T ) at a distance of 1000 km from the 

00 O 

ridge and for plate velocities of 10, 5 and 2 
cm/yr. For a given plate velocity, the lower 
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branch is for both non~Newtonran rheologies 
of Table 1 as well as for the Newtonian I 
rheology and the upper branch is for the New- 
tonian II rheology. These results are for 

3 

V*=ll cm /mole, except of course for the 
Newtonian II case. 

Figure 5. Effective viscosity t \{ 21 ) and vertical velocity v 

profiles for a dry olivine rheology and for the same 
T^, u^, V* and age values as in Figure 2. The 
top part of the figure is for = 2 cm/yr and 
the lower part is for = .10 cm/yr. Vertical 
velocity is positive upward. For a given age, 
e.g. 10 My or 150 My, there are two solutions, 

3 

one for V* = 11 cm /mole, the other for V* = 

3 

30 cm /mole, specifically indicated for the 
10 My case at the top of the figure. For the 
other- Tsets of solutions, the respective posi- 
tions of the 11 and 30 cm /mole curves are simi- 
lar . The shaded areas only approximately 
represent the range of solutions for intermedi- 
ate V* values since solution curves for differ- 
^ent V* values intersect. 

Figure Effective viscosity n (27) and vertical velocity v 
profiles for a dry olivine rheology with T =1200 
and = 2 cm/yr (top) or 10 cm/yr (bottom) . 

Vertical velocity is positive upward. The 
format is sxmilar to that of Figure 5 and the 
solutions are the n and v profiles corresponding 
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to the T and u profiles of Figure 3. The shaded 
regions suggest the range of solutions for values 

3 ' 

of V* between 11 and 30 cm /mole. 

Figure 7. Comparison of temperature T and horizontal 
velocity u profiles for wet and dry olivine 
rheologies for ages of 10 and 150 My and for 
= 1200 u^ = 10 cm/yr, V* — 30 cm^/mole. 

The temperature profile at 150 My for simple 
boundary layer cooling is also shown for refer- 
ence (dotted curve) . The dotted curve for 
10 My is practically identical to the temper- 
ature profile for the wet olivine case (solid 
curve at 10 My) . 

Figune 8. Temperature T and horizontal velocity u profiles 
for a Newzonian rheology with parameter values 
given in Table 1, T^ = 1200 °C and u^ = 10 cm/yr. 
The niombers on the figure indicate ages of 10 My 
and 150 My. For each age 3 solutions are shown. 
Solid lines correspond to the Newtonian I case 
and two values of V*. The lower T and u pro- 

3 

files are for V* = 11 cm /mole., the upper ones 
for V* = 30 cm^/mole. The shaded area between 
these curves indicates the range of' solutions 
for intermediate values of V*. The dashed 


lines are solutions for the Newtonian II case. 



Figure 9 . 


Viscosity profiles for a Newtonian rheology. 
Three solutions are shown for each of the ages 
10 My and 150 My. They correspond to the cases 
illustrated in Figure 8 . For the Newtonian I 
rheology (solid lines) and a particular age, 
the curve with the lower minimum viscosity 

3 

corresponds to V* = 11 cm /mole, the other one 

3 

to V* = 30 cm /mole. Again the area between 
these two curves is indicative of solutions for 
intejrmediate V* values. The dashed curves are 
the viscosity profiles for the Newtonian II 
rheology. The inflection of the 10 My curve 
below the region of minimum viscosity is due 
to the large shear heating bxamp on the geotherm 
for that case. Remember however that the rapid 
increase of n below the asthenosphere is arti- 

F 

f icially imposed by the' particular boundary 
condition we choose for great depth, i.e. T = 

T . 

CO 

Figure 10. Shear stress as a function of sub-asthenospheric 
temperature (T^ - T^) for different sets of 
rheological parameter values and u^ = 10 cm/yr, 

3 

X = 1000 km. a - dry olivine, V* = 11 cm /mole; 

3 

b “ wet olivine, V* = 11 cm /mole; c - Newtonian 

3 

II; d - Newtonian I, V* = 11 cm /mole. 
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Figure 11. Lithospheric thickness (depth at which n/u^ = 

. i' 

0-95) as a function of (age) Solid curves - 
dry non-Newtonian, dashed curves - Newtonian I, 
long-short dashed curve - Newtonian II. The 
Newtonian curves are for u^ = 10 cm/yri The 

3 

upper dashed curve has V* = 11 cm /mole, the 

3 

lower dashed curve has V* = 30 cm /mole. For 

3 

the upper two solid curves V* = 11 cm /mole, 
and for the lower two solid curves V* = 30 

3 

cm /mole. The value of u^ distinguishes the 
solid curves with the same value of V*; the 
thicker lithosphere is for u^ = 2 cm/yr, 
the thinner one is for u^ = 10 cm/yr. 

Pigjure 12. Surface heat flow and shear stress as functions 

• of age. Heat flow den?eases while shear stress increases 
wixh. age. Results are for non-Newtonian rheolo- 
gY/ V* = 30 cm^/mole and T^ = 1200 °C. a - dry 
olivine; b - wet olivine. The curves in the 
upp^ part of the figure are for u^ = 2 cm/yr, 
those in the lower part are for = 10 cm/yr. 

The dotted line is the surface heat flow vs. 
age for simple boundary layer ^ cooling. It is 
a shraight line with slope 

Surface heat flux and shear stress as functions 


Figure 13. 


of age. T^ = 1200 solid curves - Newtonian I, 
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dashed curves - Newtonian II. 30 and 11 are 

3 

values of V* in cm /mole. The dotted straight 
line is the simple boundary layer cooling, 
surface heat flow vs. age dependence; its 
slope is 

Figure 14. Heat flow at the ocean floor as a function of 
age for a wet olivine rheology with T =1600 °C 

3 

and V*=ll cm /mole (solid line) , and for a dry 
olivine rheology with T^=1200 °c and V*=30 cm^/mole 

3 

(long-dashed line) and V*=ll cm /mole (short-dashed 
line) . In all cases the plate velocity is 10 cm/yr. 
Heat flow data are from Sclater et al. (1976) 

(shaded rectangles) and Sclater and Francheteau 
(1970) (filled circles) . 

Figure 15. Bathymetry of the ocean floor as a function 

j, 

of (age) ^ for a dry olivine rheology with 
T^ = 1200 °C, u^ = 10 cm/yr and V* = 11 and 

3 

3-0 cm /mole. The shaded regions represent 
the bathymetry solutions for intermediate 
values of V*. In the upper and lower parts 
.or the figure a is 3x10 ^ K ^ and 4x10 ^ K 
respectively. The dashed lines are the ocean 
floor depths for simple boundary layer cooling. 

The data points are from Parsons and Sclater 
(1975) . 

Figure 16. Seismic velociides (v^ - shear, v^ - compressional) 
vs. depth for 10 My and 150 My ages. For a 
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gxven age, the profile with the lower minimum 
value of velocity is for = 1600 while 
the other is for T = 1200 ^C. The shaded area 

CO 

be-tween two such solutions indicates the lo- 
cations of velocity profiles for intermediate 
values. These profiles are based on the 
elastic properties of pure dry olivine 
(Kumazawa and Anderson, 1969) and are calculated 
from the geotherms for dry olivine, u^ = 

2 cm/yr (solid lines) . The dashed lines are 
profiles for = 1200 = iO cm/yr and 

V* = 30 cm /mole. For 10 My. the dashed lines 

coincide with the T = 1200 °C solid lines. 

/ 

Figure 17 . Comparison of geotherms based on simple 

boundary layer cooling for 10 My and 150 My 
and = L200 °C and 1600 °C, with various 

i 

solidus curves, a - dr.^ solidus? b - solidus 
in presence of hydrated minerals; c - wet 
solidus. These three solidi are from Kushiro 
et al. (1968) . The dashed solidus is that of 
I to and Kennedy (1967) for dry peridotite. 
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